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& ' Abstract 

S: 

£Nj ' In this note we evaluate Wilson coefficients for "deep inelastic scattering" (DIS) in M = 4 

SYM theory at NLO in perturbation theory, using as a probe an R-symmetry conserved 

j— ,, current. They exhibit uniform transcendent ality and coincide with the piece of highest 

transcendentality in the corresponding QCD Wilson coefficients. We extract from the 

JL ■ QCD result a NNLO prediction for the M = 4 SYM Wilson coefficient, and comment 

q_) . on the features of its Regge limit asymptotics. 

& 

CN ! 1 Discussion 

>. 
(N; 

Among the many ways through which the outstanding simplicity of the N = 4 super Yang-Mills 
r- — , \ (SYM) theory reveals itself, the pattern of transcendentality exhibited by many of the observables 

j." ■ computable in a close analytical form has the merit of setting potentially a quite direct link to 

QCD. In particular, the maximum transcendentality principle [1, 2] (MTP) is the conjecture - 
inspired by special properties [3] for the maximally super symmetric generalization of BFKL and 
evolution equations - that in the anomalous dimensions of leading twist operators only terms of 
highest transcendentality arise, which can be picked up by the 'most complicated' terms of the 
^ corresponding QCD results [4, 5] with the appropriate color factor prescription Ca = Cf = N c 

and TfUf = 2N C . Here a transcendentality weight n is given to each Riemann £ value C, n = C( n ), 
with a similar tallying for the harmonic sum Sfi(j), and the principle states that the anomalous 
dimension j(j) at n loops is a linear combination of harmonic sums of transcendentality 2n — 1. 
Several signals of consistency [6, 7, 8, 9, 10, 11] have lead to assign a predictive power to the 
MTP, which (combined with other QCD-related properties [59]) has long been the computational 
strategy for extracting multi-loop anomalous dimensions of twist operators from algebraic Bethe 
equations [12, 13, 14, 15]. 

Remarkably, patterns of uniform leading transcendentality - in this case the degree of transcen- 
dentality is 2n for n-loop results - appear also in M = 4 SYM scattering amplitudes [11, 16] 4 
even at the subleading-color (non-planar) level [18], in light-like Wilson loops [19] 5 as well as in 
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4 See also [17] for an example of how transcendentality imposes strict limitations on the admissible functional 
forms of the soft anomalous dimension. 
5 See also [20]. 



form factors [21, 22, 23] and correlation functions [24, 25] 6 . 

Significantly, the as yet heuristic nature of the observations and the presence of some exception to 
the rule of "direct extraction" from QCD results 7 make the detection of such generalized MTP 
an always necessary, preliminary check to set its predictive power in the context at hand. 

We have here evaluated the simplest non-trivial (NLO) contribution to "deep inelastic scat- 
tering" Wilson coefficients for J\f = 4 SYM, using as external hard probe the current conserved 
under the internal S77(4) R-symmetry of the theory. Although in a conformal theory bound- 
state "hadrons" do not form, one can imagine to redefine the unphysical asymptotic states 8 via 
quark/gluon distributions, governed by the J\f = 4 SYM DGLAP equations [1], in which all emerg- 
ing collinear divergencies are factorized out 9 . Looking at the final result (2.7)-(2.9), of immediate 
evidence are a) its extreme simplicity (all rational terms in z cancel and only logarithms of z and 
1 — 2 survive), b) its uniform transcendentality-degree two, evident in its Mellin transform (2.11)- 
(2.13), and c) the fact that it coincides with the highest transcendental part of the QCD result 
(2.4)-(2.5), provided the appropriate set of splitting functions reflecting the field content of the 
theories is taken into account and with the color factor prescription mentioned above. These are 
not a priori obvious observations. At least in the case of anomalous dimensions for leading(two)- 
twist operators - which govern the IR-divergent 0(l/e) 10 part of the inclusive DIS cross-section 
- the maximal transcendent ality property seems to be deeply connected to the maximal super- 
symmetry of the gauge theory (MTP is broken already in the J\f = 2 case ). While in this sense 
the 1/2 BPS [32] current (2.6) is the best candidate to gain properties a) and b) for the 0(l/e) 
contribution 12 it is non-trivial that they extend to the IR-finite part. It is also non-obvious that 
the maximal transcendentality observed for J\f = 4 SYM form factors [22] would be maintained 
in the cross-sections here analyzed. Interestingly, using the same basis (in Mellin space) which 
diagonalizes the twist-two anomalous dimensions of M = 4 SYM [1], the result simplifies further 
and can be written - in total analogy with the anomalous dimension case - in terms of a uni- 
versal function C un i{j) expressed as combination of harmonic sums and appearing with shifted 
arguments in the first two entries of the Wilson coefficient vector (3.4). Even more interesting 
are the vanishing of the third vector component, and a shift in the first entry which cancels the 
singular behavior in the Regge limit. Elaborating on these features and extrapolating them to the 
next order in perturbation theory we propose - extracting it from the QCD result [33] - the DIS 
Wilson coefficient for M = 4 SYM at NNLO, verifying its self-consistency with the analysis of its 
Regge (J — > 1) and quasi-elastic (j — > oo) asymptotics . 

Since it usually requires two non-trivial orders to reliably identify the general structure of a 
perturbative computation, it would be important to extend the calculations to NNLO and confirm 
the predictions (3.10) and (3.22), and with them the properties described above. Furthermore, 
one could proceed generalizing this computation to crossing-related processes such as the super- 
symmetric generalization of the e + e~ annihilation and Drell-Yan lepton-pair production [34]. 



See also the field-theory related structure observed for superstring amplitudes [26]. 
7 In the case of scattering amplitudes, see the mismatch [27], mentioned in [18], between the terms with highest 
transcendentality in the leading-color two-loop QCD amplitude of [28] and the corresponding M — 4 SYM amplitude. 

8 See also [29]. 

9 In an operative sense, this logic is not different from what one experiences in QCD, where the true, physical, 
asymptotic states are not reachable via perturbation theory and are actually defined via partonic distribution 
functions. 

10 We use the FDH dimensional reduction scheme of [30, 31]. 

"Private communication by L. Lipatov. 

12 It is not difficult to check that using the R-symmetry singlet, non conserved, current jp-y^ip would add to the 
formulas obtained, together with a UV divergent part due to bubble diagrams, also terms of transcendentality one 
and zero. 



This note proceeds with the presentation of the main result in J\f = 4 SYM, its comparison to 
QCD (Section 2) and the extraction from QCD of an NNLO prediction (Section 3). Appendices A, 
B and C recall, respectively, the LO expressions for splitting functions appearing in our analysis, 
a convenient basis for expressing our result and our prediction and basic definitions of harmonic 
sums. 



2 DIS Wilson coefficients in M = 4 SYM at NLO 

In the QCD deep inelastic scattering analysis, the Wilson coefficients C q , C g are the short-distance 
functions which appear in the IR-safe part of the structure functions, where they multiply respec- 
tively the (renormalized) quark and gluon distributions. At order 0(a s ) and in MS factorization 
scheme, the structure function F^ reads 

F 2 (x,Q 2 )=^^/ 1 ^,g,^)F g (^,,) + x^^/ 1 ^ 5 g,^)F 9 g,,) (2.1) 

q,q Jx " q,q Jx " 

2, 



where x is the Bjorken variable x = Q /(2p ■ q) and 



F q (9., z) =S(l-z) + ^ (p®>(z) log 9. + C^ D {z)) + 0(a% 



(2.2) 

Mf >*) =S( p - 0)(2)log S + ^ QCDw ) + ° ( ^- (2 - 3) 



Here \i is the factorization scale [35], the LO non-singlet splitting functions Pqq ,Pqg are reported 
in Appendix A and C q , C g are given by [36, 37] 



C^{z) = C F [-2p qq (z)(\og( T ^) + l)+^ + ^z-(9 + 4C 2 )5(l-z)], (2.4) 
C 9 QCD (z) = T R \ - 2 p qg (z) ( log (-^)+ 4) +6}. (2.5) 



Above, p qq ,p qg are the polinomials (A.4)-(A.5) usually introduced in literature and representing 
the highest transcendental part of the splitting functions. As a result of factorization, both 
splitting functions and Wilson coefficients C q , C g can be evaluated perturbatively squaring and 
integrating over the phase space the relevant form factors. At NLO, the processes contributing to 
C q are of real gluon emission 7* + q — > q + g and virtual gluon correction 7* -f q -> q, while the 
initial gluon process g + 7* — > q + q contributes to C g (see for example [38] ) . 

To mimic the electromagnetic interaction in the M = 4 SYM theory case [39, 40] we will 
substitute the electromagnetic current with the SU{4)r conserved current 

fJ = ^T I ^il>--ci)T I (-i^ + il^, 1 = 1,..., 15, (2.6) 

where a summation over the SU(N C ) and SU(4)r indices is understood. Notice that, although 
with the choice above we are intentionally not selecting a U(l) subgroup of the internal SU(4)r 
symmetry so that all the fermions and scalars of the theory are "charged" on equal foot 13 , the 
absence of truly non-abelian structures in the actual computation makes the effect of the sum 



13 In [39], the U(l) subgroup generated by the diagonal generator £3 = diag(l, —1,0,0) was selected, under which 
only two of the Weyl fermions and two of the complex scalars are charged. 



over all the 15 generators T 1 in the vertex with the u SU(4)ji photons" practically equivalent to 
the consideration of 15 U(l) photons. 

Together with the presence of many more processes due to scalar contributions both in virtual 
and real diagrams, the main difference with respect to the QCD setting stays in the necessity of 
assigning appropriate R-symmetry (flavor) factors to the squared diagrams and interference terms. 
This is due to the presence, in the Yukawa-type vertices of J\f = 4 SYM lagrangian, of sigma- 
matrices connecting the vectorial and spinorial representations of SO(6)r ~ SU(4:)fi, under which 
the M = 4 SYM scalar and fermionic fields, respectively, rotate. In performing the calculation, 
we have used Feynman rules derived from the M = 4 SYM lagrangian in the notation of [41] and 
the FDH regularization scheme of [30, 31]. 

The "DIS" Wilson coefficients of M = 4 SYM are then given at NLO (order a = ^=) by the 
following formulas, extremely simple and symmetric in scalar and fermionic contributions, 



Cx(z) 
C g (z) 



i(0) 



2C A \P^\z) log (— ) + 2C 2 C (/) <5(1 - z) 



2C A 



Pf (z) log 



,(0) 
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2C A [pW(z) log (— ) + 2( 2 C {V) 6(1 - z) 



where P is defined as the following linear combination of splitting functions 



(0), 



Pi 0) (z) = P£! (z)C {f) + Pj (z)C (v) a = g^,X 



(2.7) 
(2.8) 
(2.9) 

(2.10) 



Above, the LO splitting functions are those reported in Appendix A, formulas (A.9)-(A.ll), Ct v \ 
and Cft\ are the quadratic Casimir invariants respectively in the vectorial representation of SO (6) 
and fundamental of SU(4) (the latter is equivalent to the spinorial representation of SO(6)), C A 
is the Casimir eigenvalue for the adjoint representation of the gauge group SU(N C ). Already 
in this form, it is not difficult to see that the first two functions C\ and C g above, with the 
appropriate truncation of the scalar sector, coincide with the highest transcendental part of their 
QCD counterpart (2.4)-(2.5) 14 . In Mellin space, formulas (2.7)-(2.9) read 



C x (j) = 2C {f) C A 2 5 2 (j) + 4 



Mi) 



C g (j) = 2C (f) C A 



CJj) 



2C , („) C A 



4 Ml 
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25 2 (j) + 3 
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SiU) 
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2Sn(j) 



2S u (j)--2 
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(2.11) 
(2.12) 
(2.13) 



Assigning to -^ transcendentality degree n, the expressions above exhibit uniform transcendental- 
ity of degree 2, which is the maximum degree expected at this order in perturbation theory 15 . In 
that formulas above are not only given in terms of harmonic sums of fixed degree 2, their uniform, 
leading transcendentality feature is sometimes referred to as "weak", as opposed to the one in 
"strong" sense which is the argument of next section. 



14 For the definition of the degree of transcendentality for a general function, see for example [42, 43] and references 
therein. 

15 It is well-known [44] that the maximum order of soft and collinear divergences of the Sudakov form factor at 
one loop is \. 



3 Change of basis and prediction at NNLO 

Experience in the calculation of anomalous dimensions [1, 2] of leading twist operators in M = 4 
SYM suggests that the n-loop anomalous dimension matrix 7 a b(j) (a, b = g,X, (j>), after diagonal- 
ization, assumes an intriguing form in terms of a single universal function J un i(j) expressed only 
through a combination of harmonic sums of constant degree 2n — 1. Explicit expressions for the 
matrices entering the diagonalization can be found in Appendix B. It is quite natural to ask how- 
Wilson coefficients would appear in this new basis. 

Let us first notice that the particular combination of splitting functions (2.10) appearing in 
formulas (2.7)-(2.9) would also appear as coefficient of log (%) in the N = 4 analogous of (2.2)- 
(2.3). In Mellin space it is therefore natural to define the vector 

TOO = (%U),7x(J),%U)), (3-1) 

with 16 7a(j) = i\aU) + §7ia0')> a = 9i^i^i an d rotate it using the same matrix (B.3) used 
in [1, 2]. Setting v g = 2y-£^ and v\ = j^j in (B.l) and (B.3), the result is 

7(i) = V-'lij) = (l ( l(j - 2),72(i),0) , (3.2) 

where 7„ n j(x) = — 4Si(a;). Turning to the Wilson coefficients it is natural to define 

C(j) = (C g (j),C x (j),C (j> (j)), (3.3) 

where the three components are given in (2.11)-(2.13) 17 . Rotating again with V~ 1 given in (B.3) 
one gets 

C(j) = V- l C{j) = fc uni (j - 2) - jp^yi, Cum(j), o) (3.4) 

where we defined 

C uni (j)=A(S 1 , 1 (j)-S 2 (j)). (3.5) 

Comparing this result with the QCD one-loop coefficient functions in Mellin space [33], one realizes 
that the maximally transcendental part of that result is given precisely by the combination (3.5) 
of harmonic sums of degree 2. It is therefore for the Wilson coefficients rotated vector that 
the proper maximum transcendentality principle works, in total analogy with [1]. For the case 
presented here, however, a couple of new interesting features appear. The main difference with 
respect to the structure of the diagonal expression (B.2) is the vanishing of the third component 
for the vectors (3.2) and (3.4). It is quite interesting that the same feature is present in both cases 
and it would be interesting to check it at NNLO. The other novel feature is the shift of , -_iy2 in 
the first component (3.4). It is worthwhile to notice that this additional term is responsible for 
the cancellation of the singularity appearing in C un i{j — 2) in the Regge limit j — > 1. 

To better understand the pole structure in the Regge limit, it is interesting to consider the 
QCD case [45, 33, 46], where the gluon Wilson coefficient Fg (1, j) has the following form 



i? CD (u) = i + E {^§r) ef } ' QCD (i) , 

n=l ^ ' 



(3.6) 



16 The coefficient | comes from the ratio -Jr^- obtained after collecting C(/j in the Mellin transform of (2.10). 

17 

From now on an overall factor C^j^Ca will be neglected. 



which, for j — >■ 1, reads 



C^ CD (j) = U f + 0((j-l) 1 ), 



CM><**>U) = \n f C A 



43 



2C2 



1 



+ o((j-i) 



CP^ C %) = f n,tf 



1234 



13C 2 +4C 3 



1 



(j ~ 1)' 



+ o((j-ir 1 ), 



(3.7) 



with nf the number of active quarks in QCD. Starting from NNLO (n = 2) above, the most 

singular terms in Cg (j) contain (j — l) 1 "™ contributions and are proportional to the factor nf. 

Direct inspection of the QCD results shows that, both for anomalous dimensions and coefficient 
functions, the contributions with a factor nf do not contain terms with maximal transcendentality, 
2n, expected at a given loop-order n. So, according to the MTP, the terms ~ rif, including the 
singular ones in (3.7), should not contribute to the M = 4 SYM case. Accepting, by analogy 
with the case of anomalous dimension matrix, that at all orders in perturbation theory there is 
one universal Wilson coefficient function obeying the MTP, it is easy to see that the only possible 
poles in its '+' component 18 - i.e. the first one of the vector in (3.4) - are combinations of terms 

with 2 < m < 2n — 1. Since such terms are absent in the QCD 



2n—m 



(j), we conclude 



l/(j - l) In and Cm/ {J ~ 1) 

case, and C„"j(i) should contain only the most complicated terms of Cg 

that no singular terms should appear in C\ (J) (see eqs. (3.21) and (3.22) below). This would 
imply that all the singularities for j > — > 1 are contained in the first component of the anomalous 
dimension matrix (or, equivalently, in the first component of vector (3.2)). 

Another argument supporting the finiteness of the coefficient function in the Regge limit comes 
from the analysis of BFKL equation for J\f = 4 SYM. The correspondence between the singularities 
for j — > 1 obtained by the BFKL approach and those appearing in the "+" component of the 
anomalous dimension (i.e. "f U ni(j — 2)) has been checked up to NLO and it is believed to hold 
at all orders in perturbation theory [1]. In other words, the "+" component of the diagonalized 
anomalous dimensions matrix in QCD and in M = 4 SYM can be reconstructed from the BFKL 
approach in the limit j — Y 1 (see [1, 47, 3]) 19 . Hence, assuming that the BFKL equation collects all 
singular terms, for j tol, of the J\f = 4 SYM structure function F un i{Q 2 /u 2 ,j), all the singularities 
should appear only in the anomalous dimension 7 U m(j — 2) and not in the Wilson coefficients. 



All these arguments suggest that the "+" component of the Wilson coefficient rotated vector 
should be finite in the Regge limit. It would be interesting to check this conjecture at NNLO by 
explicit computation. 

Following these considerations we can guess the form of C un i(j) at NNLO from the correspond- 
ing QCD result [33]. The structure of the result in Mellin space, in analogy with (2.1), should 
be 



r ,'91 ■ 

r uni \ 2 ' J 






+ a' 



l+'a 
^(i) + f^(i)+7S(i)C«(j))l<. 



Q : 



(3.8) 



+ ^o-))V 



+ Oia 



pM 



18 In analogy with the anomalous dimension case, it is the first component of the rotated Wilson coefficient vector, 
in which Cuni(j — 2) appears, which exhibits singularities in the Regge limit. 

19 Several authors speak about the reconstruction of gluon anomalous dimension, but this is not true beyond LO. 
What is reconstructed is the first component of the diagonalized matrix (see the recent discussion in [48]). 



where 20 

2 C uni = S l ~ S 2, (3-9) 

\ Cuni = Sf-8 Sl (s 2 + S_ 2 ) + 4 5x (s_ 3 + 2 3_ 2)1 ) + 5 2 2 + 4 S 2 S„ 2 + 6 S 2 _ 2 

+10S_ 4 + 2 54 + 853,i-85- 3j i-45_ 2:2 . (3.10) 

The definition of the nested harmonic sums S± a> ±t, i ± C) ...(j) is given in Appendix C. The extraction 
of C„ ni from the QCD result [46] is rather complicated, and here we will report only the expression 
for j — y oo. In this limit and up to three loops, one finds 21 

C«(j) = 21og 2 f -2C2 + Oir 1 ), (3.11) 



C^(j) = 21og 4 ^-8C 2 log 2 ^- 16C 3 log4- + ^C4 + 0(j- 1 ), (3.12) 

Jo Jo Jo 2 

C2(i) = t lo g 6 -- 12 ^log 4 --32C3log 3 - + H9C4log 2 ^- 
3 Jo Jo Jo Jo 

-(352C5- 2 fC 2 C3)log^ + ^C3 2 -^C6 + 0(j- 1 ), (3.13) 

with jo = e~ jE . 

Let us comment on the agreement of expressions (3.11)-(3.13) with the expected results coming 
from exponentiation of the large j logarithms 22 . From general arguments [51, 52, 53] it is known 
that in DIS, for a general massless gauge theory, large logarithms for j — » oo are controlled by 
the cusp anomalous dimension r cusp (d) and an additional anomalous dimension TDis(a). In the 
case of J\f = 4 SYM, with vanishing /3-function, including finite contributions for j — > oo along 
the lines of [54] one gets 

C um (j)~g(a)e G (^), (3.14) 

where g(a) collects all the finite contributions of C U ni(j) for j — > oo. For the simple case of N = 4 
SYM the function G ( d, -I- | can be written as 

G (a, L j = \ r cusp (a) log 2 I- + T DIS (a) log 4- . (3.15) 

V Jo/ 2 jo Jo 

The anomalous dimensions and the function g(a) clearly admit a perturbative expansion in powers 
of d. From our result (3.11)-(3.13) we can read out the first terms of those expansions, getting 

g(a) = 1 - 2 d C2 + f d 2 C 4 + d 3 (^Cs 2 - ^Ce) + 0(d 4 ) , (3.16) 

r cU s P (d) = 4d - 8d 2 C 2 + 88d 3 C 4 + 0(d 4 ) , (3.17) 

iWa) = -16d 2 C 3 + 32d 3 C|c 2 C3 + HCsl + 0(a 4 ) . (3.18) 



20 Note that, as in Ref. [47, 49], our normalization of y(j) contains an extra factor —1/2 with respect to standard 
literature (see [60]) and differs by a sign in comparison with [4]. 

21 We are grateful to G. Korchemsky for pointing out a misprint in the three-loop result in the first version of the 
paper. 

22 See also [50] for a prediction on the large j behavior of the universal structure constants of TV = 4 SYM, which 
can be compared - invoking the OPE analysis of DIS - with our predictions (3.11)-(3.13). 



It is noteworthy that the function r cusp (a) coincides up to three loops with the well-known expan- 
sion of the cusp anomalous dimension of TV = 4 SYM 23 [9, 11]. This provides a good compatibility 
check for our result. 

Let us now consider the limit of C^(j) for j — ► — l 24 . To this purpose it is useful to recast 
the coefficient functions (3.9)-(3.10) in the form 25 

Cffiitf) = *(Si,i-S 2 ), (3.19) 

c uniti) = 48^1,1,1,1 " 56 ( 5 2,i,i + Si,2,i + 51,1,2) + 56^3,1 + 48S 2 , 2 + 24 (s 2 ,_ 2 + S_ 2 ,_ 2 

-16 (Si-2,1 + 25i,i,_ 2 ) + 8 5_ 3 ,i + 40 (,§1,3 + 5i,- 3 ) - 12 5_ 4 - 28 S 4 , (3.20) 

The asymptotic behavior of the expressions above can be extracted from Appendix A of [2], 
putting r = — 1. The assumption of finiteness as j — > 1 for the "+" component of the vector C 
leads to the following conjecture for its explicit form at NNLO 

eft!) = CiUO - 2) - y^yj, (3.21) 

cfO) = cfflo-q-^ + ^^fc-^Lfc. (3.22) 

It would be interesting to verify this conjecture by explicit computation. 
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A Splitting functions in QCD and in M = 4 SYM 

The leading order contributions to the DGLAP splitting functions in QCD read 

P${x) = C F {2p qq (x) + 3S(l-x)}, (A.l) 

P${x) = 2n fPqg (x), P$(x) = 2C FPgq (x), (A.2) 

P$(x) = ±C APgg (x) + 11Ca+ 6 AnfTR 5(l-x) (A.3) 



23 Our convention for the cusp anomalous dimension follows the one used in [9]. The function 7fc(a) defined in 
[11] is given by Jk(a) = 2r cU s P (a) with a — 2a. 

24 This is equivalent to study the limit j — V 1 of C^iij — 2), which is the function appearing as the first component 
of (3.4). 

25 As j — > — 1, nested harmonic sums S± a ,±b,...{J) have a less singular behaviour than standard sums S± a (j) [2]. 
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Their corresponding, Mellin-transformed, anomalous dimensions are 
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The splitting functions for twist-two operators in J\f = 4 SYM are given by [1] 
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P^(z) = C A 4(l-2z + 2z 2 ); 
P^{z) = C A A; 



P$(z) = C A 3z; 
P$(z) = C A 6(z-z 2 ); 
Pi%) = C A: ' 
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Their corresponding anomalous dimensions are 
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(A.9) 
(A.10) 
(A.ll) 

(A.12) 

(A.13) 
(A.14) 

(A.15) 
(A.16) 



B Diagonalization of anomalous dimensions in j\f = 4 SYM 



In [1] it was seen that the matrix J a b(j) (a, b = g,X, cfr) is diagonalized by 

/ v Q 



V 



-2(j-lK (,-+1)^2) ^ 



V V 9 l(3+l)Vq 



3+i v <t> 



Namely 



V'^V 



lumij - 2) 

7uni(j) 

7«™(j' + 2) 



(B.l) 



(B.2) 



with 26 

(l + 1)(l+2) (& + 3 K l 4^(2i + 3)«ji 3 u i^ aJ (2i + 3)t,J 1 

2(47 2 -l)(27+3) J i + l W 9 D (i+1)0'+2) V 9 6 (j + l)(j+2) V q I 

7 \ (2j - I)*" 1 -4(2j - l)^ 1 3(2j - l)^ 1 

The explicit expression of j un i up to NNLO is given in [2] . 

C Harmonic sums 

The basic definition of standard and nested harmonic sums with general indices are 

S ±*W = E ~^> S ±"^,~(J) = £ ^S ±6 , ±c ,...(m), (C.l) 

m=l m=l 

where we omit the sign "+". The nested harmonic sums appearing in (3.20) are defined as 

S-a,b,C-(J) = ("I)' S- a ,b,c,..Xj) + S- aAc ,...(oo) (l - (~iy) , (C.2) 

S-a-b,c,-U) = S -a,-b,c,...(j) + U - (-l) J j5'_ biCi ...(oo) 5_ a (oo) - S-aO') , 

^,-6,c,..(i) = (-1)" 5 a ,_ 6)C ,...0-) + (l - (-l) j ) (Sa,-6,c,...(oo) - 5_ 6>Cl ...(oo)[5«(oo) - S„(j)])i 

Sa,b-c,-(j) = (-1) J S at b-c„.(j) + [I - (-1) J J (S ,6 ,- c ,...(oo) - S_ C) ...(oo) S a ,b(oo) ~ S a ,b(j) 

- \Sb-c,-(oc) - S b (oo)S- C} ...(oo)) S a (oo) - S a (j) J. (C.3) 

Expressions (C.2)- (C.3) are defined for integer values of the arguments (see [55, 56]), but can be 
easily analytically continued to real and complex j via the methods of Refs. [1, 57, 56, 58]. 
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